Landau-gauge propagators in Yang-Mills theories at (3 = 0: 



massive solution versus conformal scaling 
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We study Landau-gauge gluon and ghost propagators in Yang-Mills theories at lattice parameter 
f3 = 0, considering relatively large lattice volumes for the case of the SU(2) gauge group in three and 
four space-time dimensions. We compare the lattice data to the so-called massive and conformal- 
scaling solutions, examining the requirements for a good description of the propagators over various 
ranges of momenta and discussing possible systematic errors. Our analysis strongly supports the 
massive solution, i.e. a finite gluon propagator and an essentially free ghost propagator in the infrared 
limit, in disagreement with Ref. [lj. Moreover, we argue that discretization effects play no role in 
the analysis of these propagators. 
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I. INTRODUCTION 

In recent years, considerable effort has been invested in 
the study of the infrared (IR) behavior of Green's func- 
tions in Yang-Mills theories. The results of these stud- 
ies, using analytic methods as well as numerical lattice 
simulations, are usually compared to the predictions of 
various confinement scenarios. Since Green's functions 
are gauge-dependent quantities, one can expect to find 
different confinement pictures when considering different 
choices of the gauge. 

Here we consider the Landau gauge and test the pre- 
dictions of the Gribov-Zwanziger and Kugo-Ojima con- 
finement scenarios at lattice parameter (3 — 0, i.e. in 
the limit of infinite lattice coupling. This study is very 
similar to that presented in Refs. 1, 2, 3]. On the other 
hand, it is not intended to be a simple duplication of that 
work since, in particular, the analysis we present here is 
not done exactly in the same way as in Ref. [lj]. Also, 
even though there is probably no real difference at the 
level of the lattice data, as we show in Section III CI and 
explain in Sections II El and IIIII we mostly do not agree 
with the data analysis and the interpretation of the re- 
sults presented in Ref. [l|. The interested reader is of 
course invited to read both papers and form her/his own 
opinion on this subject. 



A. Confinement scenarios in Landau gauge 

In Landau gauge, the Gribov-Zwanziger scenario 4. 
IE IS relates confinement of quarks to a ghost propa- 
gator G(p) enhanced in the IR limit when compared to 
the tree-level behavior 1/p 2 as a function of the momen- 
tum p. Indeed, in this scenario, the enhancement of G(p) 
should account for the long-range mechanism responsible 
for confinement. An enhanced ghost propagator is also 
obtained by the Kugo-Ojima color-confinement criterion 



[1, [j|. Consequently, in both cases one expects to find 
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G(p) = oo, which is referred to as ghost domi- 



At the same time, the gluons should be confined due 
to the violation of reflection positivity [12] • This implies 
that the gluon propagator in position space D{x) should 
be negative for a range of values of the space-time sepa- 
ration x. Since the gluon propagator at zero momentum 
D(p = 0) is proportional to J d d xD{x), it is clear that 
negative values for D(x) will tend to reduce the value of 
D(p = 0), leading to a suppression of D(p) at small mo- 
menta |4|. Il2l. [L3J. Similarly, in the Kugo-Ojima scenario 
one can show that the perturbative massless pole of the 
gluon propagator probably disappears as a consequence 
of the confining criterion Q. Thus, an IR-suppressed 
gluon propagator can also be accommodated in the Kugo- 
Ojima confinement scenario [Io| . Clearly, maximal vio- 
lation of reflection positivity is obtained if D(p = 0) = 0. 



Even though these two scenarios predict similar behav- 
iors for the gluon and ghost propagators in Landau gauge, 
one should recall that the line of thinking in the two cases 
is quite different. Indeed, in the Gribov-Zwanziger sce- 
nario confinement is due to the properties of the config- 
urations belonging to the boundary of the so-called first 
Gribov region [J, [1, 0] , which should be the relevant con- 
figurations in Landau gauge. On the contrary, in the 
Kugo-Ojima scenario 8, 9], confinement is obtained if 
one can define unbroken global color charges. Moreover, 
in the latter case, Gribov copies do not seem to play 
any role in the confining mechanism. Nevertheless, if one 
uses the Kugo-Ojima confinement criterion as a boundary 
condition [141 ] then the partition function is equivalent to 
the one obtained using the Gribov-Zwanziger approach, 
i.e. by restricting the functional integration to the first 
Gribov region. 
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B. Massive and conformal-scaling solutions 

Yang-Mills theories can be studied non-perturbatively 
using the field equations of motion of the theory, i.e. the 
sets of coupled Dyson-Schwinger equations (DSE) [i~5| . 
Usually, the solution of these equations depends on the 
considered truncation scheme and on the approximations 
employed. 

Recently, it has been shown [HI [13] that there exist two 
possible consistent solutions of the DSE in 4d Landau 
gauge. One solution, usually called conformal scaling, 
gives an IR-enhanced ghost propagator G(p) ~ p~ 2 ( 1+K °) 
and an IR-fmite gluon propagator D(p) ~ p 2 ( 2K z _1 ) with 
infrared exponents k g = n z = k € [1/2,3/4]. Note 
that .D(O) = for k > 1/2. The second solution, called 
massive or decoupling solution, is characterized by a tree- 
level-like ghost propagator at small momenta G(p) ~ p~ 2 
and by a finite nonzero gluon propagator D(0) > at 
zero momentum (i.e., kg = 0, kz = 1/2). 

These two solutions have indeed been obtained by sev- 
eral groups [H, [H, S HH H3, [M HI, using different 
truncation schemes. Let us recall that the conformal so- 
lution has also been obtained (in Landau gauge) in 2d 
and 3d Yang-Mills theories (ill l25| with the exponent re 
approximately given by 0.2(d — 1), where d = 2,3 and 
4 is the space-time dimension. Possible explanations of 
the origin of these two solutions have been discussed in 
[261 . I271 28]. For a description of confinement based on the 
conformal solution see 0, HI] . On the contrary, in [3(| 
the authors relate massive (respectively massless) gluons 
to color confinement (respectively deconfinement). Color 
confinement in the presence of massive gluons has also 
been related to the condensation of vortices [3l|. It is 
interesting to note that a criterion for quark confinement 
obtained in Ref. [HJ is satisfied both by the scaling and 
by the massive solution. 

The massive solution can also be obtained in 3d and 4d 
[33l [H, HH by using the Gribov-Zwanziger approach. In 
this case this solution appears when a suitable mass term 
is added to the action (while preserving its renormaliz- 
ability). In particular, the massive behavior is related 
to the condensation of a mass-dimension d — 2 operator. 
By setting to zero the value of this condensate, one gets 
back the conformal solution. It is interesting that the 
same approach cannot be extended to the 2d case [3||. 
Finally, the massive solution is also found by considering 
a mapping (in the IR limit) of the Yang-Mills action onto 
the Xcj) 4 theory (37|. 

C. Lattice studies 

Lattice simulations allow a true first-principles study 
of the IR sector of QCD, with no uncontrolled approx- 
imations. However, when studying the IR behavior of 
Green's functions in a given gauge, care is needed in or- 
der to control possible finite-size, Gribov-copy and dis- 
cretization effects. 



Recent numerical results for Landau gauge at very 
large lattice volumes [H, [H, 53] indicate a finite gluon 
propagator D(p) at zero momentum and a tree-level- like 
IR ghost propagator G(p) in 3d and in 4d. In particular, 
a flat ghost dressing function — or, equivalentl y, a null 
exponent kq — has been seen first in [38|, H^, |41| ] and 
more recently in [42|, |43|, [44j . A gluon propagator with 
an IR exponent kz — 0.5 has also been obtained using a 
tadpole-improved anisotropic lattice action |45l| . On the 
contrary, in the 2d case, one finds [12, IH, |47|] D(0) = 
and G(p) ~ p- 2 ( 1+K G) with Kq ~ Kz between 0.1 and 0.2. 
Note that this implies an IR exponent k in relatively good 
agreement with the prediction of the conformal solution 
[25j ]. i.e. k — 0.2. Thus, lattice simulations suggest a 
massive solution in 3d and in 4d, while the 2d case seems 
consistent with the conformal solution. It is important 
to note, however, that in the three cases one finds a clear 
violation of reflection positivity [H, |49|, HH for D(x) 
at x ~ 1/to. Also, lattice data confirm that, in the limit 
of large lattice volumes V, the measure of the functional 
integration gets concentrated on the boundary of the first 
Gribov region, in agreement with the Gribov-Zwanziger 
confinement scenario [1, [H, 0] ■ Indeed, in this limit the 
smallest nonzero eigenvalue of the Faddeev-Popov oper- 
ator goes to zero [3, [52|, [H, HH in 2d, 3d and 4d Landau 
gauge. 

The extrapolation of the gluon- and ghost-propagator 
data to the infinite-volume limit has been recently im- 
proved by considering rigorous upper and lower bounds 
[42l l47l [55l ] for -D(O) and G(p s ), where p s is the small- 
est nonzero momentum on the lattice. These bounds are 
valid at each lattice volume V and must be extrapolated 
to infinite volume, just as for the propagators. However, 
the bounds are written in terms of quantities that are 
easier to compute, better behaved or more intuitive than 
the propagators themselves. This allows a more precise 
extrapolation and may provide a clearer interpretation of 
the behavior of the propagators in terms of statistical av- 
erages. We note that similar bounds can also be written 
for D(p) and G(p) at general lattice momentum p, and 
for various gauge-fixing conditions. Thus, the bounds 
can be used to check the necessary conditions for the IR- 
enhancement of G(p) and for the IR-suppression of D(Q), 
clarifying when a Gribov-Zwanziger-like confinement sce- 
nario can be considered for a given gauge [56| . 

Gribov-copy effects on gluon and ghost propagators 
have been extensively studied on the lattice (43l . l57l . l58l . 
[59l . [oOj . Recent results [HI, using an extended gauge- 
fixing procedure, suggest that the restriction of the con- 
figuration space to the so-called fundamental modular re- 
gion A 5] produces a slightly more IR-suppressed gluon 
propagator. At the same time, all studies agree that 
this restriction makes the ghost propagator less singular. 
This result has a very simple explanation [57] if one re- 
calls that the fundamental modular region A is a subset 
of the first Gribov region f2 and that the smallest nonzero 
eigenvalue \ m in of the Faddeev-Popov matrix M. is larger 
for configurations belonging to A than for configurations 
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in f2 |52j. Then, since G(p) ~ M.^ 1 and \ m in goes to 
zero in the infinite-volume limit, it is natural that G(j>) 
should be smaller when a restriction to the fundamental 
modular region is implemented, i.e. for any finite lattice 
volume we must have G\(p) < Gn(p), where the sub- 
scripts indicate the region considered for the functional 
integration. Of course, after taking the infinite-volume 
limit one would still find G\(p) < Gn(p). 

Finally, discretization effects are important for the 
breaking of rotational symmetry as well as for the possi- 
ble different discretizations of the gluon field and of the 
gauge-fixing condition. In order to reduce effects due to 
the breaking of rotational symmetry, three different ap- 
proaches have been considered: to cut out the momenta 
characterized by large effects [6l[ (the so-called cylin- 
drical cut) , to improve the lattice definition of the mo- 
menta 62] or to include (hypercubic) corrections into the 
momentum-dependence of the Green's functions [63l.[&^|. 
As for the discretization of the gluon field and of the lat- 
tice Landau gauge condition, several different definitions 
can be considered [H, [U M, III [H, IM [Z^. These 
studies have usually found that different discretization 
procedures lead to gluon pro pagators that differ only by 
a multiplicative constant |66l . l69l . Fn| . The situation may 
be different for the (3 = case, discussed below. 



D. The case (3 = 

The case (3 = corresponds to the (unphysical) strong- 
coupling limit of lattice gauge theory. However, when 
studying the IR behavior of Green's functions in a given 
gauge, it has some interesting advantages compared to 
the usual simulations in the so-called scaling region. In- 
deed, at (3 = the partition function determines that the 
gluon configuration should be completely random, just as 
a spin model at infinite temperature. Thus, correlation 
functions are probing only the effects due to the gauge- 
fixing condition. Moreover, since at (3 = the lattice 
spacing is infinite (see Section III Al below) , any lattice 
volume considered is also infinite and we can hope for 
an easier study of the infinite- volume limit of the theory. 
Finally, one should also recall that the inequalities ob- 
tained in [jjl, [H| for the Landau-gauge gluon field and 
gluon propagator are valid for any value of /?, including 
(3 = 0. 

Early numerical studies of the SU(2) gluon pr opag ator 
at small (3, including (3 = 0, were presented in |57lFr2l|. 
showing the first numerical evidence of a gluon propa- 
gator suppressed in the IR limit. In particular, it was 
seen that, for [3 = 0, the gluon propagator is decreasing 
(roughly monotonically) as the momentum p decreases. 
At small positive (3, the gluon propagator also decreases 
for p below a certain turn-over value pto, with p to de- 
pending on the value of (3. For this study the largest 
lattice volume considered was V = 30 4 . We note again 
that no such behavior is observed in the scaling region 
on symmetric 4d lattices. 



E. Recent results at (3 = 

A recent and extensive study at (3 = 0, for the SU(2) 
case and considering lattice volumes up to 56 , has been 
reported in [H, [3, In this case, the authors justify the 
consideration of the case (3 = because it should corre- 
spond to the formal limit of Aqcd — > oo. This, in turn, 
would allow a study of the behavior of Green's functions 
for a range of momenta n/L <C p -C Aqcd, which is a 
necessary condition for the observation of infrared behav- 
ior on the lattice. 

Quoting the Conclusions of Ref. [U, the authors find 
that: 

• the propagator's dressing functions show conformal 
scaling behavior for large lattice momenta, a 2 q 2 ^> 

i; 

• finite-size effects are small; 

• the combined gluon and ghost data are consistent 
with an IR exponent k = 0.57(3); 

• both propagators show massive behavior at small 
momenta, i.e. a 2 q 2 <C 1; 

• this massive behavior depends strongly on the lat- 
tice discretization used for the gluon field and for 
the gauge fixing; 

• while it is possible that this ambiguity disappears 
in the continuum limit, it is definitely present at 
commonly used values of the lattice couplings in 
SU(2). 

We will comment on the analysis of the data in Section 
III CI Here we would like to stress that we do not agree 
with the limit Aqcd oo as a motivation for this type 
of study. First of all, the authors of [l[ do not quote a 
value for the lattice spacing at (3 = or for Aqcd- Thus, 
it is not clear how the comparison is really made. Of 
course, fixing the lattice spacing in the strong-coupling 
regime is not a simple issue (we will comment again on 
this issue in Section III A[) . Indeed, since we are far away 
from the continuum limit, the use of an experimental 
input is not really justified. Nevertheless, if for example 
we fix a from the physical value of the string tension a, 
then the strong-coupling expansion [in the SU(2) case] 
gives [H| a 2 a = -log(-/3/4) and in the limit (3 — > 
we find a — > oo. In this case all masses, evaluated 
on the lattice at (3 = 0, will have a null physical value 
if expressed in physical units and a comparison to the 
continuum physical value of Aqcd would be essentially 
meaningless, since even the mass of a very heavy hadron 
would be infinitely small compared to Aqcd- One could 
also try to set a finite value for the lattice spacing, as 
done for example in Ref. [73[ using the results presented 
in [74| . Then, depending on the mass used as a physical 
input, one finds a lattice spacing ranging from 0.455/m to 
1.4/7n. This corresponds to an inverse lattice spacing 1/a 
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between 141 and 433MeV. Then, assuming for example 
a value for Aqcd of about 225MeV — which corresponds 
to Ajgg — we find that Aqcd would have a value between 
0.62 and 1.92 in lattice units at /3 = 0. Thus, in this case, 
the condition p <C Aqc D is not satisfied by almost all the 
momenta considered in this work and in Refs. [2, 0, ||- 
Moreover, if one selects as significant for the IR limit 
only data verifying this inequality, then the data region 
with scaling behavior (described in the first item above) 
should be discarded. 

As for the ambiguity in the results related to different 
discretizations, it seems to us that it is not really present 
if one considers all data already available [98( . Indeed, in 
Fig. 11 of Ref. [l|, the difference between the standard 
and the modified Landau gauge for the running coupling 
constant when (3 is 2.3 is about 0.84/0.73 « 1.15 at the 
smallest momentum (corresponding top 2 « 0.048 GeV 2 ), 
for which the effect seems to be larger. This is quite a 
large difference but, as we will see below, it is proba- 
bly mostly due to the ghost sector. Indeed, the ghost 
dressing function enters quadratically into the definition 
of the running coupling, i.e. this effect is artificially en- 
larged by considering the running coupling instead of 
the propagators. This difference decreases at larger mo- 
menta, being about 12.5% at p 2 w 0.060 GeV 2 , 10% at 
p 2 w 0.222 GeV 2 and 4.5% at p 2 w 0.518 GeV 2 . (Here 
we used a — 0.83814 GeV^ 1 in order to convert lat- 
tice momenta to physical units.) Now .if one looks at 
the data reported in Figure 1 of Ref. \m for the case 
P = 2.5, there are "hardly any differences" between the 
two different discretizations, as stressed by the authors 
themselves. To be more precise, the difference in the 
gluon field is clearly within error bars, while for the ghost 
dressing function the ratio between the two results seems 
close to 2.06/2.02 rj 1.025, again considering the smallest 
nonzero momentum (corresponding to p 2 « 0.206 GeV 2 ). 
This implies a difference of about 5% for the running cou- 
pling (to be compared to a 10% difference at a similar 
value of p 2 for (3 — 2.3). For the next point, correspond- 
ing to p 2 m 0.413 GeV 2 , this difference is about 2.4%. 
Thus, the discretization effects observed in Ref. [l| seem 
to disappear in the continuum limit, being already re- 
duced by a factor 2 when going from (3 — 2.3 to p = 2.5. 
Of course, a volume V = 32 4 at p = 2.5 is rather small. 
On the other hand, data for the gluon and the ghost 
propagators atpsj 0.4 — 0.6 GeV, usually do not suffer 
from large finite-size effects 71]. One should also recall 
that large discretization effects in the ghost sector when 
using the so-called modified Landau gauge can be easily 
explained: indeed, in this case the discretized Faddeev- 
Popov matrix has an extra term, which is quadratic in the 
gluon field [l], Iz3|. This term, of course, is not present in 
the continuum expression for the Landau Faddeev-Popov 
term and is also not present in the standard lattice Lan- 
dau gauge. Thus, it is very plausible that the large dis- 
cretization effects observed in [l| be due to this peculiar 
characteristic of the modified Landau gauge. This could 
also explain why in this case discretization effects cannot 



be accounted for by a global multiplicative factor as in 
[IE HI, l7l| . By considering all this, we do not find the 
discretization effects studied in [l[ to be of significant rel- 
evance in the analysis of the IR behavior of Green's func- 
tions in Landau gauge. In any case, we are now studying 
these discretization effects in more detail [75| . 

Finally, we would like to stress that while simulations 
at null P may be interesting for a qualitative analysis 
[76| . they should clearly not be taken too seriously at 
the quantitative level. Thus, a value of k close to the 
preferred value of the so-called scaling solution — see 
the third item above — should at most be considered as 
a peculiar coincidence, not as a physically relevant result. 



II. NUMERICAL SIMULATIONS AND 
RESULTS 

Here we present data for the gluon and the ghost prop- 
agators in Landau gauge at P = both in 3d and in 
4d. The 3d data, with lattice volumes up to 100 3 , have 
been presented in [77]]. The analysis presented here and 
the 4d data, at V = 64 4 , are new. We consider the 
SU(2) case. Let us recall that recent numerical studies 
[H IzE Izl S H have verified that IR Landau-gauge 
gluon and ghost propagators are rather similar for the 
SU(2) and the SU(3) gauge groups, as expected from 
DSE studies. Thus, the analysis presented here is likely 
valid for the SU(3) case as well. 

In the 3d case we considered seven lattice volumes, 
namely V = 10 3 , 20 3 , 30 3 , 40 3 , 60 3 , 80 3 , 100 3 , with (re- 
spectively) 1000, 900, 773, 700, 1240, 344 and 364 con- 
figurations. In the 4d case we produced data only for 
the lattice volume V — 64 4 , with a total of 567 con- 
figurations. The gauge fixing has been done using the 
stochastic-overrelaxation method [13]. For the ghost 
propagator, due to technical reasons related to the com- 
puters used for the simulations, we employ the point- 
source method described in [54], [83]]. Let us note that 
the use of the point-source method usually increases the 
statistical noise in the evaluation of the ghost propagator 
as compared to the plane- wave source [54| . However, in 
the 4d case with more than 550 configurations we have 
that more than 70% of the data have a relative error 
smaller than 3% and 94.5% of the data have a relative 
error smaller than 5%. So, the introduced fluctuations 
are clearly compensated by our increased statistics. 

For our study of the gluon propagator in the 3d case we 
considered data for three different types of momenta, i.e. 
momenta with components (0, 0, q), (0, q, q) and (q, q, q), 
plus possible permutations of the Lorentz index. This 
implies that data corresponding to the momenta (0, 0, q) 
and (0,g, q) have a statistics three times as large as the 
data corresponding to the momenta (q, q, q). In the ghost 
case we evaluate the propagator only for the momenta 
(0, 0, q) and (0, q, q) (plus permutations of the Lorentz in- 
dex). Similarly, in 4d we consider momenta with compo- 
nents (0,0, 0,q), (0,0, q,q), {0,q,q,q) and {q,q,q,q). For 
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the ghost (respectively gluon) propagator we fully (re- 
spectively partially) applied permutation of the Lorentz 
index. 



A. Lattice spacing at f3 — 

As shown above, it is not simple to fix the lattice spac- 
ing in the strong-coupling regime. Indeed, if one uses an 
experimental input the result varies from a finite large 
value to an infinite value. This uncertainty is of course 
a manifestation of the possible discretization effects en- 
countered at (3 = 0. 

Nevertheless, we believe that for the study of corre- 
lation functions in Landau gauge at /3 = one should 
consider the lattice spacing as infinite. This result does 
not use an experimental input but it is based on the 
properties of the gauge-fixing algorithms used for fixing 
the random configurations to Landau gauge. Since all 
the "dynamics" in the present study is produced by the 
gauge-fixing procedure, it seems to us more reliable to 
use the gauge-fixing process as an input for fixing the 
lattice spacing. To this end we notice that if a — oo then 
all lattice volumes correspond to infinite lattice size (in 
physical units) and all (non-zero) correlation lengths are 
also infinite (again in physical units). This means that 
at P = we can work at constant physics — i.e. fixed 
correlation length in physical units — by simply chang- 
ing the lattice volume [85| without tuning the value of 
/3, as one has to do when using finite nonzero values of 
(3. (The same is true at /3 = oo, where all lattice sides 
and correlation lengths are zero in physical units.) In- 
deed, the evaluation of critical slowing-down of standard 
Landau gauge-fixing alg orithms, which requires working 
at constant physics [82j . can be done also at = [8^ |. 
yielding the expected critical exponents. 

Of course, since we assume a = oo, in this work all 
data will be presented in lattice units. 



B. Breaking of rotational invariance 



In the continuum, the gluon propagator D(p) and the 
ghost propagator G(p) are just functions of p 2 = J^p^- 
On the lattice, the momentum components p M are usually 
given by p M = 2 sin(7rn At /A r M ), where is an integer in 
the interval 1 — . . . , and A M (here supposed 

even) represents the number of sites in the /j, direction. 
However, as explained in the Introduction, due to the 
breaking of rotational symmetry, one expects the gluon 
and ghost propagators on the lattice to depend also on 
hypercubic corrections [62], |6J, H3| , with the leading term 
proportional to pM = YluP%- I n this Section we try to 
quantify these corrections and see if they can introduce 
systematic effects on the analysis of the data. Note that, 
to this end, we do not need to consider the /3-dependence 
of these corrections, as described in Section II El above 




FIG. 1: The gluon dressing function p 2 D(p) as a function of 
the (unimproved) lattice momenta p 2 for the lattice volume 
V — 100 3 . Symbols +, A and O represent data correspond- 
ing (respectively) to momenta (0,0, q), (0,q,q) and (q,q,q). 
Recall that for these types of momenta the largest value of p 2 
is respectively equal to 4, 8 and 12. 



for the discretization errors of the gluon field, but only 
estimate the corrections at the considered value of (3. 

In order to get rid of these corrections one should ex- 
trapolate the numerical data to the limit pM — > 0. This 
can be done when data are available for momenta charac- 
terized by the same value of p 2 but with different values 
of p", i.e. they belong to the same 0(4) orbit but to 
different H(4) orbits (see [(33| for description and nota- 
tion). With the choice of momenta considered in the 3d 
case here (see the third paragraph of Section |TT] above) 
the extrapolation pM — > cannot be easily done, since 
we essentially have only one #(4) orbit for each 0(4) 
orbit. The situation does not improve if one uses a dif- 
ferent discretization for the lattice momenta [63|, e.g. 
p p = 2im fJi /N fJi . 

As said in the Introduction, rotational symmetry can 
also be (partially) restored by considering [62| an im- 
proved lattice momentum with magnitude 



Pi 



V 



1 + s p [4] jp 



(1) 



where s is a numerical coefficient. (When considering 
physical units, s has mass dimension —2, i.e. it is propor- 
tional to a 2 .) Note that, with the momenta considered 
here in the 3d case and for a fixed momentum p 2 ^ 0, the 
correction pM jp 2 is maximal, and equal to p 2 , if the mo- 
mentum p has only one component different from zero. 
On the contrary, this correction is minimal, and equal 
to p 2 1 d (where d is the space-time dimension), if all the 
components of p are equal. Also note that the cylindrical 
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cut [6l| selects data for which this correction is relatively 
small, i.e. data close to the diagonal direction (q,q,q). 

In order to estimate the value of the coefficient s in Eq. 
(Q]) , we consider 3d gluon data corresponding to a certain 
type of momenta, e.g. (0, q, q). Then, for a given value of 
s, we find a spline describing these data as a function of 
the improved lattice momentum p 2 . Finally, we use this 
spline as a fitting function for the data corresponding to 
the other two types of momenta — (0,0, q) and (q,q,q) 
in this example — again considered as a function of the 
improved lattice momentum p?, with the same value of 
s. Let us recall that a spline is usually quite unstable 
outside the set of data that it describes. Thus, this fit 
has been done only for data corresponding to a value of 
p 2 that is inside the range of momenta used for the evalu- 
ation of the spline. Also, the zero-momentum datum has 
not been used for this analysis. Clearly, the chosen value 
of s, for each lattice volume, is the one that minimizes 
the x 2 /dof value for the fit. We find that the coefficient 
s obtained in this way depends on the type of momenta 
used to evaluate the spline. This is mainly related to 
the fact that the ranges described by the three different 
types of momenta are quite different. However, in most 
cases we found a value for s of the order of w 0.01 — 0.02 
and in all cases we found values smaller than the per- 
turbative value s = 1/12 w 0.0833 [iH, which is usually 
a good guess for data in the scaling region (see for ex- 
ample the ghost propagator in maximally Abelian gauge 
[871]). That the effects of breaking rotational invariance 
are small is also confirmed (see Fig. [TJ by the plot of 
the gluon dressing function. Note that these effects are 
usually more visible when considering p 2 D(p) instead of 
D(p) [88j. The situation is similar for the ghost propaga- 
tor in the 3d case and for the gluon and ghost propagators 
in the 4d case. 

One should note that having small rotational- 
symmetry-breaking effects at (3 — is perhaps a surpris- 
ing result, since one might expect strong discretization 
errors given that a is infinite. This is probably due to 
the fact that at = all data are effectively in the deep 
IR region and in this limit violation of rotational sym- 
metry is usually quite small. 



C. Plots and fits 

In order to test for behavior according to the confor- 
mal or to the decoupling solutions — and in analogy 
with Refs. [H, [1, HJ — we tried a fit of the numerical 
data using for the gluon propagator the fitting function 
6 + c(p 2 ) 2Kz+1 ~ d / 2 , where d is the number of space-time 
dimensions, and using for the ghost propagator the fit- 
ting function c(p 2 ) _KG_1 . For the gluon propagator we 
considered the cases 6 = and 6^0, corresponding 
to Eqs. (18a) and (18b) of Ref. Q. The function used 
for the ghost corresponds to Eq. (19a) of the same ref- 
erence. However, note that we usually fit the ghost 
propagator and not the ghost dressing function as in 



TABLE I: Parameter fit for the gluon propagator D(p) as 
a function of the (unimproved) lattice momentum p 2 using 
the fitting function b + c(p 2 ) 2KZ ~ 0,5 and data points in the 
range p 2 £ [0, 1.5]. We do a separate fit for each of the three 
types of momenta (0, 0, q), (0, q, q) and (q, q, q), here indicated 
(respectively) as 1,2 and 3. For each fit we also report the 
number of degrees of freedom (dof) of the fit and the value of 
X 2 /dof. Notice that we considered only the lattice volumes 
V > 40 3 in order to have enough data points for the fit. 



V 


momenta 


b C KZ 


dof X z /dof 


40 3 


1 


0.319(8) 0.299(9) 0.51(2) 


6 


2.38 


40 3 


2 


0.28(3) 0.34(3) 0.47(4) 


2 


3.64 


40 3 


3 


0.33(5) 0.30(5) 0.5(1) 


1 


3.22 


60 3 


1 


0.308(6) 0.307(6) 0.50(1) 


10 


2.40 


60 3 


2 


0.29(1) 0.33(1) 0.48(2) 


5 


1.90 


60 3 


3 


0.23(6) 0.37(6) 0.42(5) 


3 


3.13 


80 3 


1 


0.314(6) 0.307(6) 0.52(1) 


14 


1.50 


80 3 


2 


0.29(1) 0.32(1) 0.48(2) 


8 


1.37 


80 3 


3 


0.309(9) 0.296(9) 0.51(2) 


6 


0.24 


100 3 


1 


0.305(5) 0.307(6) 0.50(1) 


18 


1.73 


100 3 


2 


0.28(1) 0.33(1) 0.46(1) 


11 


1.61 


100 3 


3 


0.30(1) 0.32(1) 0.49(2) 


8 


0.97 



Ref. (we will comment again on this point in Sec- 
tion III C 51 below) . Also recall that the scaling solution 
[HI gives D(p) ~ ( p 2)2«+i-d/2 and G ^ _ ( p 2 )-K-\ 

with k « 0.3976 in 3d and k « 0.5953 in 4d. The fits 
have also been done separately for the different types of 
momenta used in our simulations. This allows again an 
estimate of systematic effects due to the breaking of ro- 
tational invariance. 



1. 3d gluon propagator 

In the gluon case we find that data at small momenta 
are best fitted with b =/= 0. Then, the exponent kz is very 
close to 0.5 in the 3d case (see Table A} if one uses for the 
fit the range p 2 G [0,p^J withp 2 n = 1.5. A slightly larger 
or smaller value for kz is obtained if one uses respectively 
p 2 m — 1 or p 2 m — 2. For momenta larger than p^ the best 
fit is obtained considering 6 = 0. In this case (see Table 
HI)) the exponent kz is very close to 0.35 and this value is 
essentially independent of the value of p 2 n . In both cases 
the effects due to the breaking of rotational invariance 
are relatively small, as expected (see discussion in Section 
III 131 above) . Finally, we tried a fit of all the data with 
6 7^ (see Table IIII[) . Clearly, in this case one finds 
for kz a kind of average between the value 0.5 found 
at small momenta and the value 0.35 obtained at large 
momenta. It is interesting to notice that a value very 
close to the conformal solution k 0.3976 is obtained for 
the largest lattice volumes and for momenta of the type 
(0, q, q) and (q, q, q). As an example we report in Fig. [3] 
the fits obtained at small and at large momenta for the 
lattice volume V = 60 3 and for data corresponding to 
the momenta (0,0, q). 
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TABLE II: Parameter fit for the gluon propagator D(p) as 
a function of the (unimproved) lattice momentum p 2 using 
the fitting function c(p 2 ) 2rez_0 ' 5 and data points in the range 
p 2 > 1.5. We do a separate fit for each of the three types of 
momenta (0, 0, q), (0, q, q) and (q, q, q), here indicated (respec- 
tively) as 1,2 and 3. For each fit we also report the number of 
degrees of freedom (dof) of the fit and the value of \ 2 /dof. 



V 


momenta 


C HZ 


dof xVdof 


10 3 


1 


0.65(2) 0.34(2) 


1 


1.43 


10 3 


2 


0.640(3) 0.349(2) 


2 


0.11 


10 3 


3 


0.65(4) 0.34(1) 


2 


2.77 


20 3 


1 


0.628(9) 0.357(6) 


4 


0.85 


20 3 


2 


0.623(6) 0.355(3) 


6 


1.25 


20 s 


3 


0.64(2) 0.344(7) 


6 


2.11 


30 3 


1 


0.6291) 0.356(8) 


7 


2.12 


30 3 


2 


0.631(7) 0.351(3) 


9 


1.46 


30 3 


3 


0.625(9) 0.351(4) 


10 


0.91 


40 3 


1 


0.637(7) 0.348(5) 


10 


1.15 


40 3 


2 


0.620(4) 0.356(2) 


13 


0.61 


40 3 


3 


0.64(1) 0.347(4) 


14 


1.58 


60 3 


1 


0.625(5) 0.355(4) 


16 


1.23 


60 3 


2 


0.620(3) 0.355(2) 


20 


0.70 


60 3 


3 


0.618(7) 0.355(3) 


22 


1.65 


80 3 


1 


0.626(7) 0.354(5) 


22 


1.41 


80 3 


2 


0.627(5) 0.352(2) 


27 


0.77 


80 3 


3 


0.624(9) 0.352(4) 


29 


1.17 


100 3 


1 


0.616(5) 0.362(4) 


28 


1.06 


100 3 


2 


0.625(4) 0.354(2) 


34 


0.87 


100 3 


3 


0.629(6) 0.349(2) 


37 


0.62 



TABLE III: Parameter fit for the gluon propagator D(p) as a 
function of the (unimproved) lattice momentum p 2 using the 
fitting function 6 + c(p 2 ) 2KZ ~ ' 5 and all data points. We do 
a separate fit for each of the three types of momenta (0, 0, q), 
(0, q, q) and (q, q, q), here indicated (respectively) as 1,2 and 3. 
For each fit we also report the number of degrees of freedom 
(dof) of the fit and the value of \ 2 /dof. 



V 


momenta 


b 


c 


Kz 


dof x'/dof 


10 3 


1 


0.49(1) 


0.16(1) 


0.55(2) 


3 


2.09 


10 3 


2 


-0.2(2) 


0.9(2) 


0.33(1) 


2 


0.17 


10 3 


3 


-1(7) 


2(7) 


0.3(1) 


2 


2.74 


20 3 


1 


0.37(1) 


0.24(2) 


0.51(2) 


8 


7.84 


20 3 


2 


0.18(5) 


0.44(6) 


0.39(2) 


7 


2.26 


20 3 


3 


-0.1(2) 


0.7(2) 


0.34(3) 


7 


2.64 


30 3 


1 


0.33(1) 


0.28(2) 


0.48(2) 


13 


8.22 


30 3 


2 


0.12(4) 


0.50(4) 


0.38(1) 


12 


2.42 


30 3 


3 


0.09(6) 


0.52(7) 


0.37(1) 


12 


1.47 


40 3 


1 


0.30(1) 


0.31(1) 


0.46(1) 


18 


5.93 


40 3 


2 


0.17(2) 


0.44(2) 


0.392(8) 


17 


2.24 


40 3 


3 


0.08(7) 


0.54(7) 


0.36(1) 


17 


2.75 


60 3 


1 


0.289(8) 


0.32(1) 


0.457(7) 


28 


5.42 


60 3 


2 


0.21(1) 


0.40(2) 


0.405(6) 


27 


3.48 


60 3 


3 


0.17(3) 


0.45(3) 


0.387(8) 


27 


2.35 


80 3 


1 


0.289(9) 


0.32(1) 


0.455(8) 


38 


3.59 


80 3 


2 


0.22(1) 


0.39(2) 


0.406(6) 


37 


2.39 


80 3 


3 


0.20(3) 


0.41(3) 


0.398(9) 


37 


1.79 


100 3 


1 


0.288(6) 0.315(7) 0.460(6) 


48 


2.45 


100 3 


2 


0.24(1) 


0.37(1) 


0.413(5) 


47 


2.35 


100 3 


3 


0.20(2) 


0.42(2) 


0.393(6) 


47 


1.40 



0.9 it 




0.5 1 1.5 2 2.5 3 3.5 4 



P 2 

FIG. 2: The gluon propagator D(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volume V = 
60 3 for momenta of the type (0, 0, q). We also present the fits 
obtained at small momenta p 2 < 1.5 and at large momenta 
p 2 > 1.5 (see Tables U and HJl. 



One can also check that the gluon propagator D(x) vi- 
olates reflection positivity. Actually, for all cases we find 
that D(x — 1) is already negative. [Of course, D(x = 0) 
is always positive.] In Fig. [3] we show D(x) as a func- 
tion of the space-time separation x for the lattice volume 
V = 100 3 when considering momenta of the type (0, 0, q). 
Following Ref. [49] , we fitted the data using a sum of two 
Stingl-like [89[ propagators f(x) = fi(x) + f 2 (x) with 
fi(x) — Ci cos (hi + Xix) e~~ XiX . We found a good de- 
scription of the data (see again Fig. [3]) by fixing c\ — 
D(x = 0) = 66.6393, 6i = 0,b 2 = tt/2,A 2 = A x /3 and 
with C2 = 20(2) and Ai = 2.5(2). Also, by comparing 
these data to results obtained in the scaling region [49| . 
one can use the value of x where D(x) starts to oscillate 
around zero as an input for fixing the lattice spacing. 
From Fig. [3] this hap pens at x ss 3, giving a rs 1/m (see 
Figure 5 in Ref. [49j). This is in quantitative agreement 
with the values obtained in [73| (see Section [IE] above). 

We also consider the volume dependence of D(p = 0) 
and the bounds introduced in (47|. As one can see from 
Fig. [U these bounds are well satisfied. Compared to the 
results obtained at finite (nonzero) /3 (see Fig. 1 in [47|). 
we find that the values of D(0) are closer to the upper 
bound than to the lower one. We have also checked that 
we can extrapolate the data for D(0) and for the upper 
and the lower bounds to a finite nonzero value as well as 
to zero (see Table HV)) . However, considering the value 
of x 2 /d°f, our data seem to prefer an extrapolation to 
D(0) ^ (see also Fig. 13}. 

Let us note that the above result is in agreement with 
what observed in Ref. !)()'. i.e. in the 3d case volumes 
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FIG. 3: The gluon propagator D(x) as a function of the space- 
time separation x for the lattice volume V = 100 3 , obtained 
by Fourier transforming data corresponding to momenta of 
the type (0,0,q). For clarity, only the interval x £ [0,8] is 
represented here. We also present the fit obtained using the 
fitting function f(x) described in the text. 

up to 100 3 are not large enough in order to have a com- 
plete control over the extrapolation to V — oo at p = 0. 
This suggests that — for the zero-momentum propaga- 
tor D(0) — simulations for (3 = and for f3 values in the 
scaling region are essentially equivalent. In other words, 
the value of -D(O) seems to be related only to the ther- 
modynamic limit (V — > oo) and is not affected by the 
value of the lattice spacing. This peculiar behavior at 
p = may be related to the fact that the gauge-fixing 
condition p ■ A(p) = does not play a role in this case. 
Also, one might argue that D(0) is simply a measure of 
the zero modes of the gluon field, whose value is of course 
strongly affected by the lattice size. The situation is dif- 
ferent for p 0. Indeed, while at /3 = the propagator 
is decreasing with p for all lattice volumes and momenta, 
for values of (3 in the scaling region one sees a decreasing 
propagator only for large enough lattice volume and for 
momenta below a certain value pto- Also, at (3 = finite- 
size effects for p > are very small or null (see Fig. [5] 
below and compare to Fig. 1 in Ref. [l|), while for j3 > 
one sees large finite-size effects for p < p to . 



2. 4d gluon propagator 

We did the same fitting analysis also in the 4d case for 
the gluon propagator D(p) at the lattice volume V = 64 4 . 
As shown in Tables El El and ELfl we find D{0) w 0.45 
and Kz ~ 0.9 at small momenta (using the fit with b =/= 0), 
while the fit for p > 1.5 (with b = 0) gives D(0) = and 
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FIG. 4: The gluon propagator at zero momentum D(p = 0), 
as well as the upper and the lower bounds introduced in [47| ]. 
as a function of the inverse lattice size 1/N. We also present 
a fit to the data using the function a + b/N c (solid line) and 
the function e/N d (dashed line) (see Table ITT)) . 

TABLE IV: Parameter fit for the gluon propagator at zero 
momentum D(p = 0), as well as for the upper and the lower 
bounds, as a function of the inverse lattice size 1/N using 
the fitting functions a + b/N c and e/N d . For each fit we also 
report the value of x 1 /dof . 



fit a b c x 1 /dof 

lower bound 0.048(3) 0.5(4) L2(4) 4~22~~ 

D{0) 0.30(2) 3(2) 1.2(3) 3.42 
upper bound 0.46(3) 4(3) 1.2(3) 3.62 

fit e d X Z /dof 

lower bound 0.12(2) 0.20(4) LTo"" 

D(0) 0.74(9) 0.20(4) 13.5 
upper bound 1.1(1) 0.20(4) 12.3 



nz ~ 0.56. (The fits corresponding to the range of small 
and large momenta are shown together with the data in 
Figs. [B] and 0) These results essentially do not change 
by cutting the data in Tables [V] and I VII at p 2 = 1 or at 
p 2 = 2. However, contrary to the 3d case, our data for 
p 2 > 1.5 can also be described by the fitting function 
b + cx 2kz ~ 1 with b 0. Indeed, in this good 
fit — with x 2 /dof ~ 0.97 and dof = 23 — is obtained 
with b = 0.4(1), c = 0.14(9) and n z = 0.67(6). Finally, 
a fit using all the momenta and b = suggests a kind 
of average value for kz, i.e. kz ~ 0.7. Thus, the fits 
suggest a value for kz close to the scaling solution result 
k « 0.5953 only if one ignores the data at small momenta 
and forces the parameter b to be zero. 

As in the 3d case, the effects due to violation of ro- 
tational symmetry are clearly small. Also, violation of 
reflection positivity is observed already at x = 1 and 
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FIG. 5: The gluon propagator D(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volumes 
V = 20 3 (symbol +), 40 3 (symbol A) and 80 3 (symbol O). 
Here all types of momenta are represented for the three lattice 
volumes considered. 



TABLE V: Parameter fit for the gluon propagator D(p) as a 
function of the (unimproved) lattice momentum p 2 using the 
fitting function b + c(p 2 ) 2Kz ~ 1 and data points in the range 
p 2 € [0, 1.5] for the lattice volume 64 4 . We do a separate fit 
for each of the four types of momenta (0,0, 0,q), (0, 0, q, q), 
(0,5,5,5) and (5,5,5,5), here indicated (respectively) as 1,2, 
3 and 4. For each fit we also report the number of degrees of 
freedom (dof) of the fit and the value of x 2 /dof. 



momenta 



kz dof x /dof 



1 


0.446(4) 0.095(5) 0.85(3) 


11 


1.33 


2 


0.460(5) 0.095(6) 0.83(4) 


6 


0.85 


3 


0.41(4) 0.14(4) 0.67(8) 


4 


1.28 


4 


0.45(2) 0.08(2) 0.9(2) 


3 


1.97 



the gluon propagator in position space D(x) is well de- 
scribed by a sum of two Stingl-like propagators f{x) = 
fi(x)+f2{x) with fi(x) = Ci cos (bi + \ { x) e~ XiX (see Fig. 
[S]) ■ Note that the values of the masses Xi are not very dif- 
ferent from the 3d case. Finally, as explained in Section 
III C ll above. one can fix the lattice spacing by comparing 
these data to results obtained in the scaling region. In 
particular, considering our data at /3 = 2.2 for V = 128 4 , 
we find that D(x) « for x > 2fm. Then, from Fig. [51 
we find a ~ 1/m as in the 3d case. 



3. 3d ghost propagator 

We now consider the data obtained for the ghost 
propagator in the 3d case. [Note that in this case we 
did not evaluate the propagator for the lattice volume 



TABLE VI: Parameter fit for the gluon propagator D(p) as 
a function of the (unimproved) lattice momentum p 2 using 
the fitting function c(p 2 ) 2rez_1 and data points in the range 
p 2 > 1.5 for the lattice volume 64 4 . We do a separate fit 
for each of the four types of momenta (0,0,0,5), (0,0,5,5), 
(0,5,5,5) an d (5,5,5,5), here indicated (respectively) as 1,2, 
3 and 4. For each fit we also report the number of degrees of 
freedom (dof) of the fit and the value of x 2 /dof. 



momenta 


C K.Z 


dof X 2 /dof 


1 


0.533(3) 0.570(3) 


17 


0.76 


2 


0.558(4) 0.562(2) 


21 


0.86 


3 


0.561(6) 0.561(3) 


23 


0.99 


4 


0.531(6) 0.566(2) 


24 


0.97 



TABLE VII: Parameter fit for the gluon propagator D(p) as 
a function of the (unimproved) lattice momentum p 2 using 
the fitting function b + c(p 2 ) 2 " z_1 and all data points for 
the lattice volume 64 4 . We do a separate fit for each of the 
four types of momenta (0,0,0,5), (0,0,5,5), (0,5,5,5) and 
(5,5,5,5), here indicated (respectively) as 1,2, 3 and 4. For 
each fit we also report the number of degrees of freedom (dof) 
of the fit and the value of x 2 /dof. 



momenta 



kz dof x /dof 



1 


0.436(4) 0.103(5) 0.76(1) 


30 


1.60 


2 


0.430(9) 0.13(1) 0.70(1) 


29 


2.20 


3 


0.39(2) 0.17(2) 0.65(1) 


29 


1.38 


4 


0.41(1) 0.13(2) 0.68(1) 


29 


1.19 




0.01 



FIG. 6: The gluon propagator D(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volume V = 
64 4 for momenta of the type (0,0,0,5) an d p 2 < 1- We also 
present the fits obtained at small momenta p 2 < 1.5 (solid 
line) and at large momenta p 2 > 1.5 (dashed line). See Ta- 
bles El and EH 
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FIG. 7: The gluon propagator -D(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volume V = 
64 4 for momenta of the type (0,0, 0, g) and p 2 > 1. We also 
present the fits obtained at small momenta p 2 < 1.5 (sohd 
line) and at large momenta p 2 > 1.5 (dashed line). See Ta- 
bles El and [VH 
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FIG. 8: The gluon propagator D(x) as a function of the space- 
time separation x for the lattice volume V = 64 4 , obtained 
by Fourier transforming data corresponding to momenta of 
the type (0,0,0, q). For clarity, only the interval x 6 [0,8] is 
represented here. We also present the fit obtained using the 
fitting function f(x) described in the text with c\ = D(x — 
0) = 36.4604, bi = 0,6 2 = tt/2, Ai = 2.15(2), c 2 = 1.7(1) and 
A 2 = 0.65(3). 
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FIG. 9: The ghost propagator G(p) as a function of the (unim- 
proved) lattice momenta p 2 for the lattice volume V = 100 3 
for momenta of the type (0, 0, q). We also present the corre- 
sponding fit reported in Table IVTTll 



V = 30 3 .] As said above, we tried a fit using the func- 
tion c(p 2 )~ KG ~ 1 . Results are reported in Table lVIIII One 
clearly sees that the infrared exponent kg decreases as 
the lattice volume increases, in agreement with (38l. 
In this case the effects due to the violation of rotational 
symmetry are more evident. Indeed, the exponent kq is 
systematically smaller for momenta along the axes. This 
is probably due to the fact that along the axes one can 
consider smaller momenta, for which the ghost propaga- 
tor is less enhanced. 

By looking at Fig. [9] it is clear that this fit, with only 
one term, can be improved. In particular, considering 
the results reported in P, [3, 0] , one can try to see how 
the exponent kq depends on the value of p 2 . To this 
end we have ordered the data points by the value of p 2 
and divided them in sets of ten data points. Then, we 
did a separate fit in each interval, again using the fitting 
function c(p 2 ) _ftG_1 . Results for the lattice volume V = 
100 3 are reported in Table HX1 One clearly sees that the 
exponent kg increases with p 2 (see also Fig. [TU]) and it is 
usually larger for the type of momenta (0, q, q). The same 
is observed for the other lattice volumes. This result is 
actually well-known. Indeed, all lattice studies of the 
ghost propagator [57], HH, [Hj have found that G(p) is 
enhanced compared to the tree- level propagator 1 jp 2 at 
intermediate momenta. (We will comment again on this 
in the Conclusions.) 
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TABLE VIII: Parameter fit for the ghost propagator G(p) as 
a function of the (unimproved) lattice momentum p 2 using the 
fitting function c(p 2 ) _KG_1 and all data. We do a separate fit 
for each of the two types of momenta (0,0, 5), (0,q,q), here 
indicated (respectively) as 1 and 2. In each case we indicate 
the smallest nonzero momentum p m in- Finally, for each fit 
we also report the number of degrees of freedom (dof) of the 
fit and the value of x 2 /dof. 



V 


momenta 


Pmin C f^G 


dof x'/dof 


10 3 


1 


0.382 3.94(3) 0.260(8) 


3 


0.39 


10 3 


2 


0.764 3.9(1) 0.35(2) 


3 


1.75 


20 3 


1 


0.098 3.77(5) 0.209(7) 


8 


4.15 


20 3 


2 


0.196 3.57(6) 0.28(1) 


8 


3.06 


40 3 


1 


0.0246 3.64(8) 0.171(7) 


18 


14.2 


40 3 


2 


0.0492 3.38(7) 0.220(8) 


18 


5.59 


60 3 


1 


0.0110 3.74(8) 0.144(5) 


28 


29.2 


60 3 


2 


0.0219 3.89(8) 0.192(8) 


28 


16.7 


80 3 


1 


0.0062 3.84(9) 0.127(6) 


38 


18.4 


80 3 


2 


0.0123 3.50(7) 0.164(6) 


38 


7.02 


100 3 


1 


0.0039 4.1(1) 0.110(5) 


48 


19.4 


100 3 


2 


0.0079 3.46(8) 0.154(5) 


48 


9.91 
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FIG. 10: The ghost propagator G(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volume V = 
100 3 for momenta of the type (0,q,q). We also present the 
fits for the first and third sets of data points reported in Table 

m 



4- 4d ghost propagator 

Finally, data for the ghost propagator G(p) in the 4d 
case show again (see Table [Xj) that the exponent kg de- 
pends on the type of momenta and systematically in- 
creases when one considers momenta closer to the di- 
agonal than to the axes. Also, the exponent kg in- 
creases with p 2 (see Table IXII and Fig. ITTj) , going from 
a very small value — close to zero — at small p 2 up 
to almost 1 for the largest momenta and for momenta 



TABLE IX: Parameter fit for the ghost propagator G(p) as 
a function of the (unimproved) lattice momentum p 2 using 
the fitting function c(p 2 ) _KG_1 and various sets of ten data 
points for the lattice volume 100 3 (we report in the table the 
smallest and the largest momenta considered for each set). 
We do a separate fit for each of the two types of momenta 
(0,0,5), (0,5,5), here indicated (respectively) as 1 and 2. For 
each fit we also report the value of \ 2 /dof [the number of 
degrees of freedom (dof) is always 8]. We do not report the 
fit for the last set because the corresponding data points have 
large statistical fluctuations. 



Pmin 


pmax 


momenta 


c 


KG 


X*/dof 


0.0039 


0.38 


1 


5.0(2) 


0.070(6) 


11.2 


0.0079 


0.76 


2 


4.4(1) 


0.102(7) 


4.40 


0.46 


1.38 


1 


4.00(2) 


0.21(1) 


0.17 


0.92 


2.76 


2 


3.98(5) 


0.29(2) 


0.19 


1.50 


2.62 


1 


3.7(1) 


0.08(5) 


0.19 


3.01 


5.24 


2 


4.2(1) 


0.33(3) 


0.04 


2.74 


3.62 


1 


6(1) 


0.6(1) 


0.31 


5.47 


7.24 


2 


3(1) 


0.3(1) 


0.28 



TABLE X: Parameter fit for the ghost propagator G(p) as 
a function of the (unimproved) lattice momentum p 2 using 
the fitting function c(p 2 )~ KG ~ 1 and all data for the lattice 
volume 64 4 . We do a separate fit for each of the four types 
of momenta (0,0,0,5), (0,0,5,5), (0,5,5,9), (5,5,5,5), here 
indicated (respectively) as 1, 2, 3 and 4. In each case we 
indicate the smallest nonzero momentum p m in- Finally, for 
each fit we also report the value of x 2 /dof [the number of 
degrees of freedom (dof) is always 30]. 



momenta 


Pmin C KQ 


x'/dof 


1 


0.0096 5.4(2) 0.075(7) 


250 


2 


0.0193 4.6(2) 0.128(9) 


238 


3 


0.0289 4.2(1) 0.16(1) 


127 


4 


0.0385 4.0(1) 0.19(1) 


55.3 



along the diagonal. This can be seen also in Fig. [l"2l 
where an effective exponent kg has been evaluated us- 
ing the relation —1 + 0.5 log [G(pi)/G(p2)]/log(p2/pi), 
where G(p\) and G(p2) are the values of the ghost prop- 
agator at two nearby momenta p\ and pi- The plot is 
done as a function of the average momentum p ave — 
(pi +p2)/2. The errors have been evaluated using the 
so-called bootstrap method with 2000 samples. For clar- 
ity, we show only results with a relative error smaller 
than 50% and with a positive value for kg (this crite- 
rion selects about 70% of the data). One clearly sees 
that the effective exponent kq is monotonically increas- 
ing. The situation is very different in the gluon sector, 
where indeed a scaling solution can be used to fit the 
data. As one can see in Fig. [T2J the effective exponent 
k z = 0.25{2 + log[D( Pl )/D(p 2 )]/log(p 1 /p 2 )} is essen- 
tially constant in this case. Moreover, at small momenta 
one finds Kz = 0.5, in agreement with the massive so- 
lution, and at large momenta the effective exponent is 
still very close to this value. Note that a constant shift 
of 0.5 is built into the definition of the gluon exponent 
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TABLE XI: Parameter fit for the ghost propagator G(p) as 
a function of the (unimproved) lattice momentum p 2 using 
the fitting function c(p 2 )~ KG ~ 1 and various sets of eight data 
points for the lattice volume 64 4 (we report in the table the 
smallest and the largest momenta considered for each set). 
We do a separate fit for each of the four types of momenta 
(0,0, 0,g), (0,0, q,q), (0,q,q,q), (q,q,q,q), here indicated (re- 
spectively) as 1, 2, 3 and 4. For each fit we also report the 
value of x 2 /dof [the number of degrees of freedom (dof) is 
always 6]. 



Pmin 


Pmax 


momenta 


c 


KG 


x'/dof 


0.0096 


0.59 


1 


6.4(2) 


0.036(7) 


101 


0.0193 


1.17 


2 


5.8(2) 


0.06(1) 


114 


0.0289 


1.76 


3 


5.4(2) 


0.09(1) 


65.4 


0.0385 


2.34 


4 


5.3(2) 


0.10(1) 


18.6 


0.73 


2 


1 


5.192(9) 0.223(5) 


0.13 


1.46 


4 


2 


5.28(3) 


0.300(6) 


0.11 


2.19 


6 


3 


5.64(5) 


0.389(7) 


0.09 


2.92 


8 


4 


5.8(2) 


0.39(2) 


0.18 


2.20 


3.41 


1 


5.78(8) 


0.38(1) 


0.06 


4.39 


6.83 


2 


5.94(9) 


0.403(9) 


0.03 


6.59 


10.2 


3 


5.8(2) 


0.39(1) 


0.04 


8.78 


13.7 


4 


7.2(7) 


0.51(4) 


0.10 


3.55 


4 


1 


5.8(2) 


0.39(2) 


0.01 


7.09 


8 


2 


5.8(4) 


0.39(3) 


0.02 


10.6 


12 


3 


19(6) 


0.9(1) 


0.24 


14.2 


16 


4 


7(2) 


0.5(1) 


0.08 



Kz- Thus, numerical results in the gluon sector should 
be rather quoted as kz — 0.5, in order to convey a clear 
indication of the precision of the results. 



5. A massive fit for the ghost propagator 

The analysis above has shown that the scaling solution 
does not describe the ghost sector. Indeed, one cannot 
find a reasonably large range of momenta where the data 
can be fitted by a power-law with a given value of the IR 
exponent kg- 

On the other hand, the ghost-propagator data 
are well described by the fitting function f{x) — 
[a — folog (p 2 + c 2 )] /p 2 , recently proposed in [24] . which 
gives a free ghost propagator in the infrared limit. The 
parameter c can be interpreted as a gluon mass [24] ]. 
Also note that this function corresponds to the small- 
momentum limit of the fitting function used in Ref. [lij 
to fit the ghost data in 3d and in 4d for values of (3 in 
the scaling region. Using the fitting function f(x) above 
we obtain, in the 3d case, the parameters a = 3.96(2), 
b = 0.92(2) and c = 0.155(6) with x 2 /dof = 0.73 us- 
ing the data p 2 < 4 (see Fig. H"4"l) . Using all the data 
one finds a = 3.95(2), b = 0.94(1) and c = 0.159(6) 
with x 2 /dof — 0.72. A similar fit in the 4d case 
gives a = 5.51(2), b = 1.45(1) and c = 0.499(7) with 
X 2 /dof = 0.65 using the data p 2 < 4 (see Fig. Us- 
ing all the data one finds a = 5.44(1), b = 1.372(8) and 
c = 0.466(5) with x 2 /dof = 0.93 (see Fig.HH). 
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FIG. 11: The ghost propagator G(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volume V = 
64 4 for momenta of the type (q,q,q,q). We also present the 
fits for the first and third sets of data points reported in Table 
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FIG. 12: The effective ghost exponent kg, evaluated using 
— 1 + 0.5 log [G(pi)/G(p2)]/log(p2/pi), as a function of the 
average momentum p a ve = (pi +f>2)/2 for the lattice volume 
V — 64 4 . Here, G(pi) and G(p2) are the values of the ghost 
propagator at two nearby momenta p\ and p2- 



For a comparison with Ref. [l[ — see Eq. (19c) and 
Fig. 9 of that reference — we have also tried a global fit 
using for the ghost dressing function the Ansatz g{x) = 
e/(p 2 + l) KG . The resulting fit gives e = 3.88(3), Z = 
0.098(8), k g = 0.271(7) with x 2 /dof = 1.63 in 3d and 
e = 6.56(7), I = 0.73(2), k g = 0.434(6) with x 2 /dof = 
1.79 in the 4d case (see Fig. [T7J). Thus, the value of 
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FIG. 13: The effective gluon exponent kz, evaluated using 
0.25 {2 + log [D(p 1 )/D{p 2 )]/log(p 1 /p 2 )}, as a function of the 
average momentum p avE = (pi + p2)/2 for the lattice volume 
V — 64 4 . Here, D(pi) and D(p2) are the values of the gluon 
propagator at two nearby momenta pi and p2. 
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FIG. 14: The ghost propagator G(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volume V = 
100 3 and momenta p 2 < 4. We also show a fit using the 
fitting function f(x) = [a — fa log (p 2 + c 2 )] /p 2 , discussed in 
the text. 

X 2 1 'dof is clearly worse when compared to the logarith- 
mic behavior considered above. In particular, from Fig. 
1171 one can see that the power-law behavior cannot de- 
scribe well the "curvature" of the data, underestimating 
the data at small and at large momenta and overshooting 
the numerical results at intermediate momenta. This is 
of course not a surprise since we have clearly shown in 
the previous subsection that a single power-law does not 
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FIG. 15: The ghost propagator G(p) as a function of the 
(unimproved) lattice momenta p 2 for the lattice volume V = 
64 4 and momenta p 2 < 4. We also show a fit using the fit- 
ting function f(x) = [a — &log (p 2 + c 2 )] /p 2 , discussed in the 
text. 



describe the ghost data, unless one selects a very small 
interval of momenta. In the 4d case, this result does 
not improve if one forces the exponent kg to be equal to 
0.562, as done in Fig. 9 of fl]. Indeed, in this case we 
find e = 8.37(6) and I = 1.23(2) with x 2 /dof = 5.7 (see 
Fig. US]). 

One should also observe that the function g(x) con- 
sidered above and in Ref. [l[ is not a truly scaling solu- 
tion h(x) — s/(p 2 ) KG . Indeed, while the latter is char- 
acterized by a constant value for the exponent kq — 
-91og[/i(a;)]/91og[p 2 ] = -[p 2 /h(x)]dh(x)/dp 2 , for g(x) 
one has an effective exponent — [p 2 / g(x)]dg(x)/dp 2 = 
K G p 2 1 ip 2 + I)- This effective exponent is monotonically 
increasing with p 2 , becoming constant and equal to kg 
only for very large momenta. Thus, g(x) is trying to 
describe the lattice data, characterized by an exponent 
kg increasing with the momentum, while "suggesting" 
a possible scaling behavior, i.e. a function with a con- 
stant exponent. In fact, as we have shown above, the 
data are much better described by the massive solution 
f(x), suggested by a recent analytic study [24j. On the 
other hand, a function such as g(x) is very poorly jus- 
tified from the theoretical point of view. In particular, 
such a fitting function for the ghost dressing function 
implies a fit e/[p 2 (p 2 + l) KG ] for the ghost propagator. 
We do not see any theoretical reason that the mass scale 
\fl should affect only "part" of the power-law behavior 
of the propagator. Here we have also tried a fit to the 
ghost-propagator data using e/[(p 2 + l) 1+KG ]: the result 
is very poor, with \ 2 /dof = 95.6 [and kg = 0.182(6)]. 

Finally, that a simple power-law is not capable of de- 
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FIG. 16: The inverse ghost dressing function l/[p' 2 G(p)] as 
a function of the (unimproved) lattice momenta p 2 for the 
lattice volume V = 64 4 . We also show a fit using the fitting 
function 1/ [a — 6 log (p 2 + c 2 )] , discussed in the text. 

scribing the lattice data should also be clear looking at 
Fig. 5 of Ref. [l|, where it is evident that the exponent 
kg depends on the momentum x considered, and from 
Fig. 9 of the same reference, where one sees that the fit 
systematically underestimates the data at intermediate 
momenta. Actually, at the end of Sec. Ill of Ref. [l[ the 
authors clearly say that a true scaling solution, i.e. their 
Eq. (19a), gives a very poor description of the ghost data 
and that their preferred fitting function is g{x), which 
only reminds one of a possible scaling behavior. This is 
probably the reason that induced the authors of [1| to 
conclude in favor of a scaling solution at large momenta 
in the ghost sector. As already stressed above, we do not 
agree with this conclusion. 



III. CONCLUSIONS 

We have studied numerically the infrared behavior of 
SU{2) Landau-gauge gluon and ghost propagators at lat- 
tice parameter (3 = 0, considering 3d lattices of volumes 
up to 100 3 and 4d lattices of volume 64 4 . By carrying out 
a careful fit analysis of the proposed behavior according 
to the scaling or to the massive solutions of DSE, we find 
that our data strongly support the massive solution, i.e. 
a finite gluon propagator and an essentially free ghost 
propagator in the infrared limit p — ► 0. Moreover, the 
gluon propagator D(x) as a function of the space-time 
separation x violates reflection positivity and it is well- 
described by a sum of two Stingl-like propagators. These 
results are in qualitative agreement with data obtained 
at finite j3 in the scaling region. 
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FIG. 17: The inverse ghost dressing function l/[p 2 G(p)] as 
a function of the (unimproved) lattice momenta p 2 for the 
lattice volume V = 64 4 . We also show a fit using the fitting 
function g(x) = (p 2 + l)^ G /e, discussed in the text. 
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FIG. 18: The inverse ghost dressing function l/[p 2 G(p)] as 
a function of the (unimproved) lattice momenta p 2 for the 
lattice volume V = 64 4 . We also show a fit using the fitting 
function g(x) = (p 2 + i)°' 562 /e, discussed in the text. 



We should stress that, in agreement with Refs. [H, [1, [H , 
a scaling solution appears in the gluon sector if one ne- 
glects the data at small momenta. As explain in Sec- 
tion II El above, we do not see any reason for excluding 
those data from the analysis. In particular, discretiza- 
tion effects at small momenta are under control and the 
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large effects observed in Ref. [l( are probably only due to 
the bad scaling properties of the modified Landau gauge. 
Moreover, the value of k clearly depends on the way the 
fits are done. In particular, a value close to the preferred 
value of the scaling solution, i.e. 0.2(d — 1) in d dimen- 
sions, is obtained only with very specific and ad-hoc fits. 
In any case, we believe that the scaling solution is clearly 
excluded by the ghost sector and we definetively do not 
agree on this point with the analysis and the conclusions 
presented in [H, @, Q. Indeed, in this case, the IR ex- 
ponent k depends on the interval considered, increasing 
essentially monotonically as the momentum increases, i.e. 
it is impossible to find a decent "window" giving a con- 
stant value for n. On the other hand, we have shown that 
the data for the ghost propagator are very well described 
by a simple functiomrecently suggested by an analytic 
study presented in [24J. This function clearly supports 
the so-called massive solution. 

As for the ongoing discussion about massive solution 
versus conformal scaling, we remark that the lattice re- 
sults may be summarized as follows. 

• In 2d Landau gauge one sees conformal scaling [42], 




• In 3d and 4d Land au g auge one finds the massive 
solution [H, m S3, 52] , 

• In 4d Coulomb gauge, the transverse gluon propa- 
gator is well described by a Gribov formula, going 
to zero at zero momentum |)3|, [94], HH 

• In the so-called A gauges, which interpolate be- 
tween the Landau gauge (A = 1) and the Coulomb 
gauge (A = 0), one clearly sees [96| that the behav- 
ior of the transverse gluon propagator gets modi- 
fied when A goes from 1 to 0, becoming closer and 
closer to the behavior obtained in Coulomb gauge 
as A becomes smaller and smaller. 

The simulations cited above are essentially all done in the 
same way, i.e. in most of the cases by ignoring effects due 
to Gribov copies, by using a standard discretization for 
the lattice action, for the gluon field and for the gauge- 
fixing condition and by using one of the standard gauge- 
fixing algorithms. Recently, to rescue the conformal so- 
lution in Landau gauge, several authors have explained 
the lattice Landau data by evoking supposedly (very) 
large effects due to Gribov copies, discretization effects 
and bias related to the use of the usual gauge-fixing al- 
gorithms. Of course, one has to verify that all possible 
sources of systematic effects are indeed under control in 
a numerical simulation. On the other hand, it seems 
unlikely to us that these effects would show up in some 
cases of the simulations described above and not in oth- 
ers. For example, why would the 2d Landau-gauge case 



be conformal and not the 3d and 4d cases if the same 
code is used in these three cases? Why should Gribov- 
copy effects be so important in 3d and 4d Landau gauge 
and not in 2d Landau gauge and, even more strikingly, 
in 4d Coulomb gauge, where one would expect stronger 
effects since the transversality condition is imposed sep- 
arately on each time slice? In our view, present lattice 
data are simply telling us that the infrared behavior of 
gluon and ghost propagators in Landau, Coulomb and A 
gauge depends on the gauge-fixing condition and on the 
dimensionality of the system (as, for example, the critical 
behavior of statistical mechanical systems). We believe 
that the present challenge is to understand why this is 
the case and that the bounds introduced by us in Ref. 
[47l | and their interpretation in terms of magnetization 
and susceptibility of the gluon field could be a key ingre- 
dient in a simple explanation of present lattice results in 
Landau gauge. 

Finally a remark about color confinement. As said in 
the Conclusion of Ref. (42J, we point out that the be- 
havior of gluon and ghost propagators at very small mo- 
menta is probably not so important for the explanation 
of confinement. Indeed, why should the behavior of a 
Green function at a few MeV, i.e. for a space separa- 
tion of about 50/m, be relevant for hadron physics, since 
the typical hadronic scale is of order of If ml Let us re- 
call that in a recent paper it has been shown that 
the appearance of a linearly-rising potential is related (in 
Landau and in Coulomb gauge) to the momentum-space 
gluon configuration A{p) for p < 1 GeV . In this region 
one can indeed observe strong non-perturbative effects in 
the gluon and in the ghost propagators in Landau gauge: 
the gluon propagator violates reflection positivity and 
the ghost propagator is enhanced when compared to the 
tree-level behavior. Thus, some important predictions of 
the Gribov-Zwanziger scenario are still verified and, at 
the same time, one can try to relate the massive solu- 
tion to the requirement of color confinement as in Refs. 
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